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ABSTRACT 

We report on the discovery of a high velocity B star, HD 69686. We estimate 
its space velocity, distance, surface temperature, gravity, and age. With these 
data, we are able to reconstruct the trajectory of the star and to trace it back 
to its birthplace. We use evolutionary tracks for single stars to estimate that 
HD 69686 was born 73 Myr ago in the outer part of our Galaxy (r ~ 12 kpc) at a 
position well below the Galactic plane {z ~ —1.8 kpc), a very unusual birthplace 
for a B star. Along the star's projected path in the sky, we also find about 12 
other stars having similar proper motions, and their photometry data suggest 
that they are located at the same distance as HD 69686 and probably have the 
same age. We speculate on the origin of this group by star formation in a high 
velocity cloud or as a Galactic merger fragment. 



Subject headings: line: profiles — stars: rotation — stars: fundamental parame- 
ters — stars: early-type — stars: individual (HD 69686) 
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Introduction 



The objects in the solar neighborhood with high pecuhar velocities can always be placed 
into one of two categories: the disk-origin or non-disk origin. The normal, young, disk stars in 
the solar neighborhood share the same circular motion with the Sun, and hav e small velocity 
residuals in a Gaussian distribution with a standard deviation of 10-16 km s~^ (iGies &: Bolton 



1986 



and references therein). The disk-origin high velocity stars are the objects whose 
peculiar velocities are far beyond the velocity dispersion of normal disk stars. They are 
often called runaway stars because the high peculiar velocities of these stars are tho ught to 
be caused by two p ossible mechanisms: 1) supernova explosion in a binary system (jZwicky 
19571 : lBlaauwlll96ll ): or 2) strong gravitational interaction in a multiple st ar system that often 
occurs in a dynamically chaotic, young, dense s tellar cluster environment (jPoveda et al.lll967l : 
Leonard &: Duncanlll988l : iGaburov et al.ll2008l ). Previous studies of runaway stars show that 
both m echanisms generate runaway objects, but an early systematic survey by lGies fc Bolton 
(119861 ) suggests that the second mechanism is dominant. On the other hand, most of the 
non- disk-origin, high velocity objects in the solar neighborhood are thought to belong to 
the old halo population. The halo stars have high peculiar velocities because they do not 
follow the circular motion of the normal disk stars but move in their own very diverse orbits. 
Because of the significant difference in the ages of high velocity stars associated with these 
two origins, the runaway, disk origin is generally assumed for young objects like high velocity 
B stars. However, in this paper, we discuss the unusual case of a nearby, high velocity B-star, 
HD 69686, and we argue that it may have been born far from the Galactic disk. 

We are currently making a spectroscopic survey of field B stars to study the evolution 
of stellar rotation. One of our targets, HD 69686, was selected because it is fairly bright 
{y ~ 7.1) and is classified as a B8 star in the SAO star catalog. In this paper, we present 
what we have learned about HD 69686. At the time of writing, the target's listing in 
the SIMBAD database includes only very basic astrometry and photometry data and no 
published references. After measuring radial velocities of our program stars, HD 69686 was 
immediately singled out by having the largest radial velocity, 148 km s~^, among our sample 
B stars. Only after this, we noticed that this star also has quite a large proper motion 
{fiacosd = —86.17 ± 0.67 ma s yr~\ (is = 7.21 ± 0.42 mas yr~^, from the newly released 



Hipparcos re-measurements by Ivan Leeuwenll2007l ). Thus, HD 69686 is a previously unknown 
high velocity B star. How old is this star now? Is it a runaway star? Where did it come 
from? Is it a binary system? These are the kind of questions that we want to answer. 

In the next section, we briefly describe the instrumental setup and how we obtained 
the final spectra of HD 69686. In §3, we present the details of the procedure we used to 
determine the key parameters, such as the space velocity, the effective temperature, gravity. 
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and age of the star. In §4, we use a modern Galactic potential model to reconstruct the full 
trajectory of HD 69686 that helps us to locate its birthplace and leads us to suggest that 
HD 69686 might not be a runaway star as we thought previously. In the final section, we 
provide additional evidence implying that HD 69686 may been born in the Galactic halo. 



2. Observation 



The spectra of HD 69686 were obtained on the 2.1-m telescope at KPNO using the 
Goldcam spectrograph (with a 3072 x 1024 CCD detector, T3KC) on 2008 November 17 
and 18. The spectrograph grating (G47, 831 lines mm~^) was used in second order with 
a CUSO4 blocking filter, and this arrangement provided spectrum coverage of about 900A 
around the central wavelength of 4400 A. The slit width was set at 1'.'3, leading to a resolving 
power R ~ 2400 (FWHM~1.83 A, measured in the comparison spectra). The integration 
time for HD 69686 was 55 s yielding spectra with S/N ~ 280 in the continuum regions. 
During the entire observing run, we took comparison (HeNeAr lamp) exposures for each of 
our program stars, including HD 69686, to ensure accurate wavelength calibration of the 
spectra. The accuracy of wavelength calibration is estimated to be around 6 km s~^ in 
velocity space by checking the wavelength fitting residuals of the comparison spectra and by 
comparing multi-night spectra of all our program stars that are not spectroscopic binaries. 
We obtained the final reduced spectra by going through the standard IRAF0 CCD image 
reduction (subtracted the bias level, divided the flat images, removed cosmic rays, and fixed 
the bad pixels/columns) and the long-slit stellar spectrum extraction procedures (traced and 
binned the spectrum, calibrated wavelength). 



3. Stellar Parameters 



The procedure to determine the physical parameters of HD 696 86 is similar to that 
used in our previous studies of cluster B stars ( Huang fc Gies 2006a| jbl). The V sini value 
was derived by fitting synthetic model profiles of He I A4471 and Mg II A4481, using realistic 
physical models of rotating stars (con sidering Roche geomet ry and gravity darkening). The 
details of this step are described in iHuang fc Gied (l2006al ). The best fit of He I A4471 
is illustrated in Figure 1. With V smi at hand, we then determined both the effective 
temperature (T^ff) an d gravity (log (7) of the star by fitting the H7 profile (see details in 
Huang &: Giesll2006bl ). The best fit of the H7 profile is shown in Figure 2. By shifting the 



http:/ /iraf.noao.cdu/ 
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best fit profiles in wavelength to match the observations we also obtained an accurate radial 
velocity of each night's spectrum. The radial velocities were transformed into the heliocentric 
frame by removing the orbital motion of the Earth (using the RVCORRECT function in 
IRAF), and these velocities are given in Table 1. We do not see any significant change in 
the radial velocity between the two nights that is larger than the errors (6 km s~^). Even if 
HD 69686 were a long period, spectroscopic binary system, any orbital velocity variations will 
be small compared to its radial velocity. Thus, we adopted the mean of the radial velocities 
in Table 1 as the systemic radial velocity of HD 69686 for later analysis (Table 2). For a 
rotating star such as HD 69686, the derived gravity (logg) represents an average of gravity 
over its visible hemisphere. It may not be a good indicator of evolutionary status of the star 
because the effective gravity on the equatorial region is lowered by centrifugal force induced 
by stellar rotation. Instead, we use the gravity at the poles of the star (log^fpoiar) to estimate 
the age of the star since the logt^poiar value is not significantly infl uenced by rotation. The 
method to estimate logf^poiar is described in Huang fc Gies ( 2006bl) and has been applied t o 
many of our studies on the evolution o f B stars (Huang fc GiesI 2008 : McSwain et al. 2008). 



Based on our model simulation results ( Huang fc Gies 2006bl ). except for those extreme cases 



where the stars spin at or very close to the breakup velocity, the estimations of log (^poiar are 
quite accurate (the statistical errors are 0.03 dex or less). 

We need one additional step to correct T^s and logg derived from the H7 flt. We 
noticed that we had to use a model with much higher temperature (16200K) to flt the 
He I A4471 line while our H7 line flt suggests that the effective temperature of HD 69686 
should be around 14500K. This discrepancy implies that HD 69686 is a helium strong star. 
When we use models with solar abundances (H:He=0. 9:0.1) to synthesize the H7 proflle 
and to flt the observed spectrum of a He-peculiar star, we get systematic errors in the 
derived parameters (Tpff and lo g g) . The systematic errors can be corrected using the results 
given by lHuang fc Gied (120081 ) who used the same technique to determine the photospheric 
parameters. The correction procedure was done in an iterative manner. We flrst determine 
the He abundance using the derived Tes and lo g q values, then corre ct Tes and log g based 
on the He abundance according to Table 3 of iHuang &: Gied (120081 ). We then determine 
the He abundance again with the updated T^s and logg values, and repeat the steps until 
convergence is reached. At the end, we obtained the flnal corrected Tgfj and logg (listed in 
Table 2) for HD 69686, with abundances typical of a He-strong star, with a number ratio of 
H:He= (0.82 ± 0.04) : (0.18 T 0.04). 

The best available parallax a nd proper motion data of HD 69686 come from the Hippar- 
cos catalogue (Ivan Leeuweru 120071 ). However, the parallax measurement, 2.35 ±0.58 mas (Icr 
range: 341~565 pes), is still too uncertain to provide us with an accurate tangential velocity 
for tracing the star's trajectory in space. We have to flnd an alternative method to obtain 
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a reliable distance of the star. Our method is described as follows: 1) by com paring the de- 



rived Tgfj and log (^poiar of HD 69686 with evolutionary tracks of stellar models (ISchaller et al. 



19921 ). we can estimate its stellar mass; 2) we then can calculate the stellar radius from the 
mass and logf^poiar! 3) from the radius and surface temperature, we can calculate the true 
luminosity of the star (we considered a distorted shape of a rotating star in this step, see 
details below); 4) with the true luminosity and the V magnitude of the star, we can derive 
the distance for assumed values of extinction and bolometric correction. The B — V color 
index from the Hipparcos/Tycho catalogu^ (ESA 1997) is — 0.15±0.01 which is basically no 
different from the zero-reddening color of a star with temperature of 14,800 K. Thus, we 
can safely se t Ay = in the last step. The bolometric correction data we used here are 
from iBalonal (119941 ). Finally, we estimated the age of the star fr om its mass and pola r grav- 
ity by comparing with the non- rotating theoretical models by Schaller et al. (Il992). The 
whole procedure was a l so rep eated with other denser model grids by iLejeune fc Schaerer 
(I2OOII ) JPemarque et al.l (120041 ) , and iMarigo et al.l (120081 ) , and we found no noticeable differ- 
ence in the derived parameters between these model grids. All the derived parameters are 
listed in Table 2. 



Heger fc Langen (120001 ): iMeynet fc Maederl (120001 ) and, more recently, lEkstrom et al. 



J2008h 'EKS models, hereafter) have predicted that, due to rotationally induced mixing 



fast rotating stars may evolve on MS tracks quite differently from non-rotating stars. Fresh 
hydrogen can be brought down to the core due to rotationally induced mixing and make 
rotating stars have a longer MS lifetime. To investigate the difference in ages between 
rotating and non-rotating stars, we calculated the ages of the rotating models with various 
spin rates at a fixed evolutionary status (at log^fpoiar = 4.04). We choose the 3 Mq model 
(the rotating model with mass closest to HD 69686) from EKS models. The result is given in 
Table 3. As expected, fast rotating stars are evolving more slowly than slow- or non-rotating 
stars with the same mass in term of surface gravity at poles. As shown in Table 3, a very 
fast rotating star with Q/Qcrit = 0.9 probably takes 20% more time than a non-rotating star 
to reach log^fpoiar = 4.04. 

It is worthwhile noting that HD 69686 is a rotating star with a moderate V sin i value 
(= 141 km s~^). Because of the centrifugal force on its surface, the projected disk of HD 69686 
on the sky is not perfectly round, and the projected angular area is expected to be slightly 
larger than tiR^^^^J D"^ {D is the distance) by a factor /, which depends on the inclination 
angle i. Because we do not know the exact value of the inclination angle, we can only estimate 
it by averaging all possible values of i. The range of i is between arcsin(V^sini/V^rit) and 



^http://archive. ast.cam.ac.uk/hipp/hipparcos. html 
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7r/2, where V^nt is the breakup velocity of the star {Vcrit = 411 km s ^ for HD 69686). We 
found the statistical mean of the acceptable inclination range is i = 60°. If we assume that 
HD 69686 (with M = AA2Mq, log^poiar = 4.04, and, therefore, i?poiar = 3.32Rq) has an 
inclination of i = 60°, then we obtain R^q = 3.51Rq and / = 1.07 from calculations based on 
simple Roche geometry. With estimates in hand for the effective temperature and gravity, 
we can compare the observed and model spectral energy distributions (SED) of the star to 
determine its angular size, and then use our estimate of its physical size to derive a second 
determination of distance. In Figure 3, we plot the optical/IR flux points that are converted 
from the B, V, I g, and JHK<^ magnitud es of the star and from the UV data from the ESRO 



TD— 1 satellite ([Thompson et al.l Il978l ). The fit between the calculated SED (Ay = 0) 
and photometry data is very good, and the normalization factor yields an equivalent, limb 
darkened angular diameter of 0.085 ± 0.004 mas. Using the projection factor / and polar 
radius given above, we arrive at a distance of 377 ± 16 pc, which is fully consistent with the 
spectroscopic distance derived earlier. 

We also apply our II7 fit method to a set of realistic models of HD 69686 to investigate 
more closely the difference between the statistically estimated logt^poiar and the true values 
of log(7poiar for modcls at various inclination angles. The model star has a mass of 4.42 Mq. 
At different inclination angles, we synthesize the H7 profiles by considering Roche geometry 
(controled by both rotational and gravitational potential on the surface) and the gravity 
darkening effect (Tiocai oc g^i^)- We adjust Tpoiar and -Rpoiar of the model stars for each 
inclination angle so that the measured values (Tmsr and logt^msr from H7 fit) always equal 
the values given in Table 2. All results are given in Table 4, where we can see that the 
range of log (Rpoiar, 4.04 ± 0.04, covers all inclination angles from i = 90° down to i = 30°. 
Below i = 30° (i.e., Vcq > 282 km s~^), the model star spins close to the critical speed 
(around i = 20°). Because of l ow probability that the star rotates with higher V^q (> 280 



km s~^) (IHuang fc Gied l2006al ) and because the corresponding inclination range is narrow 
(20° < i < 30°), the chance of HD 69686 spinning at a low inclination angle and having a 
higher log (Rpoiar beyond the range of 4.04 ± 0.04 is low. Even if HD 69686 does spin at a low 
inclination angle, its age would be older than that its polar gravity implies because it rotates 
much faster than at a high inclination angle (see Table 3). In this numerical simulation, we 
have also checked the projected area of the models on sky. We found that the projected 
area of the models at all investigated inclination angles are basically same (different within 
±1%). As the inclination angle goes lower, the / factor gets bigger. Meanwhile the model 
star rotates faster and faster (increasing V^q). In order to keep the measured logg constant, 
we have to lower the polar radius, -Rpoiar- Thus, the net effect is that the projected area on 
the sky, nR^^-y^^/D'^ multiphed by /, stays almost unchanged. In other words, it is safe to 
derive the distance using the projected area at one inclination angle. 
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It should be noted that, as we were deriving the surface log (7 of HD 69686 from H7 
profile fit, we did not consider possible systematic errors in the synthesized profiles. Our 
synthesis code took th e specific intensity profile data generated by SYNSPEC4^if|, which is 
using the VCS theory (IVidal et al.lll973l ). If the systematic error causes an underestimation 
of surface gravity, log (7, by 0.1 dex or more, HD 69686 might be younger than we estimate 
above, and could be born on the Galactic plane (see section 4). 



Many previous studies (ITobin fc Kilkennv 



Ramspeck et al.l l2001al : iLynn et al. 



200 



I 



1981 



Martin 



Saffer et all 119971 : iMagee et all I2OOI 



2OO4J ) have pointed out that some old 



evolved objects in the halo, lying on a post-horizontal branch (PHB) or a post-asymptotic 
giant branch (PAGE), could have similar atmospheric parameters (temperature and gravity) 
as young massive OB stars of population I. Could HD 69686 be an old PHB or PAGB halo 
object with the same T^q and logg values listed in Table 2? One key difference between the 
old blue halo objects and the young OB disk stars is that the former have much lower mass 
(< IMq). If HD 69686 had a mass lower than IMq, then for its observed gravity, it would 
have a small radius and a distance of less than 180 pc, which would yield a parallax more 
than 5cr higher than the Hipparcos parallax (2.35 ± 0.58 mas). Thus, we are confident that 
HD 69686 is a young, massive, main sequence B star instead of an old blue halo star with 
similar atmospheric parameters. 



4. Tracing Backward to the Birthplace 

Now we have very accurate information that we list in Table 2 about the position and 
the space velocity of HD 69686 (resolved into three orthogonal components relative to the 
Sun: the radial velocity Vr and the tangential velocity in the directions of increasing right 
ascension Vra and declination Vr,ec for epoch J2000). With help of a properly selected 
Galactic potential model, we should be able to follow the motion of HD 69686 backward 
to its birthplace. Among the many Galacti c models available in liter ature, we decided to 



choose the axisymmetric model proposed by lDehnen &: Binneyl (Il998bl ) (their model 2, DB2 



here after). Their model mass distribution was constructed to meet the constraints imposed 
by all the high quality kinematical data that were available at that time. The gravitational 
acceleration is related to the potential by l^{R,z) = — V$(-R, 2). We obtained both the 
potential and the acceleration data of the model in a format of 3D grids in space (0.1 kpc 
per step in R and 0.01 kpc per step in z), and we made a bilinear-interpolation in this grid 
to determine the acceleration at each time step. 



•^httpi/Znova. astro. umd.edu/Synspec43/synspec. html 



- 8 - 



The star's true space velocity in the Galactic potential field is its velocity relative to 
the Local Standard of Rest (LSR, hereafter), {U,V,W), plus the Galactic circular veloc ity of 
LSR, (0, Vcirc{Ro) ,0) . Here Rq (= 8 kpc) is the Sun's distance to the Galactic center (IReid 
19931 ). and Kirc(-Ro) =217 km s~^ is the circular velocity of DB2 at Rq (very close to the 
lAU recommended value, 220 km s~^). 

In order to trace the star's motion back in time, we use the reversed space velocity of 
HD 69686, -{U,V + Kirc(i?o),M^) as the initial velocity. The initial location of HD 69686 
in the Galaxy is determined from its location rela tive to the Sun and the Galactic location 
of the Sun, 8 kpc away frorn the Galactic center (lReidlll993l ) and 8 pc above the Galactic 
plane (IHolmberg et al.l 119971 ). We calculated lf{R,z) first, and then moved the star a very 
small step assuming ~^{R, z) stays constant during this small move. We updated the star's 
velocity and acceleration at the end of each step. Repeating this calculation until reaching 
the estimated age of the star (actually reaching to the zero age of the star because we move 
back in time), the integration follows the star back to its birthplace. The numerical errors 
of this integration depend on the size of each small step. We found that a step of 2 kyr 
is a good choice which results in only a very small difference (< 1 pc) from the integrated 
positions using a smaller step (1 kyr) over 80 Myr. 

The reconstructed trajectory of HD 69686 is plotted in Figure 4 as a thick solid line. 
The dotted lines are those tracks that result by changing by 1 o" the space velocity component 
VaA (with the dominant velocity error, la = 8 km s~^). Because of the very large peculiar 
velocity of HD 69686, it has a non-circular orbital motion (the circular orbit of the Sun in 
the X — Y plane is plotted as a dashed line in Fig. 4 for comparison). Currently, HD 69686 is 
quickly passing through the solar neighborhood and moving away from the Galactic center 
(U = —191 km s^^). When we first derived U, we suspected that this star might come from 
the very inner part of the Galaxy. After the full trajectory was reconstructed, we realized 
that HD 69686 has traveled a very long distance across the Galaxy before reaching its current 
location (see the top panel of Fig. 4). The birthplace of HD 69686 is located in the outer 
part of the Galaxy (about 12 kpc from the Galactic center, 4 kpc further out than the Sun). 
At its birth (73 Myr ago), HD 69686 was more than 10 kpc away from the Sun (compared 
to its current distance of 380 pc). About 26.6 Myr ago, HD 69686 reached its minimum 
distance from the Galactic center 4.4 kpc). 



Martini (120061 ) studied more than 60 high Galactic latitude B stars and he suggested 



that al l the young popula tion I B stars in his sample a re likely ejected from the Galactic 
plane. iMagee et al.l (12001! ) and iRamspeck et al.l (l2001bl ) found very similar results in their 
investigations of high-latitude blue objects. If the cause of the high peculiar velocity of 
HD 69686 is similar to other runaway OB stars, then its birthplace is expected to be in or 
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very close to the Galactic plane because star formation activities are very rare at places far 
from the plane. However, our calculated trajectory suggests that HD 69686 is an exception. 
In the lower panel of Figure 4, HD 69686's projected trajectory on the X — Z plane clearly 
shows that it was located well below the Galactic plane iz ~ —1.8 kpc) 73 Myr ago. Today 
HD 69686 (only 150 pc above the Galactic plane) is on its way plunging into the Galactic 
plane {W = —57 km s^^). Moving backward in time, it reached the highest position of its 
trajectory (860 pc above the Galactic plane) 19.5 Myr ago, and crossed the plane around 
37 Myr ago, thanks to the stronger gravitational force in the z direction in the inner part 
of the Galaxy. However, further back in time, HD 69686 was below the Galactic plane and 
further away from the Galactic center, and, therefore, it experienced a weaker 2;-force. The 
z-force was so weak that the star's trajectory began 73 Myr ago at a position well below the 
plane {z ~ —1.8 kpc) and with a velocity approaching the plane {W = +22 km s~^). The 
errors in the velocity components and the age estimation of HD 69686 are not large enough 
to make the range of its possible birthplace extend to the Galactic plane. 



5. Discussion 

The reconstructed trajectory of HD 69686 implies that HD 69686 might be born in a 
location with a particularly low density (1.8 kpc below the Galactic plane and 12 kpc away 
from the Galactic center), where star formation activity is expected to be very rare. Of 
course, we assume here that HD 69686 is a single-star system. Though our two consecutive 
night observations do not display any variation in radial velocity above the noise level, we can 
not completely rule out that HD 69686 is a binary system with a longer period (> 10 days). 
The mass-transfer events that often occur in close binary systems could totally invalidate 
our age estimation based on single-star evolutionary tracks. If HD 69686 is in a binary 
system and was ejected from the Galactic plane due to the supernova (SN) explosion of its 
companion, the explosion occured most likely about 37 Myr ago when HD 69686 was crossing 
the Galactic plane. The coordinates of the crossing position are {x,y,z) = (2919,-4477,0) 
pc, and the space velocity of HD 69686 at that moment is (V^, Vy, V^) = (—326.2, —24.4, 84) 
km s~^ while the Galactic circular velocity at that location is (-181.8,-118.5,0) km s~^ 
(using DB2). Then we can calcula te HD 69686's p e culiar velocity (i.e. ejection velocity) at 



that location, which is 192 km s ^. iNelemans et al.l (119991 ) investigated the relation between 



the ejection velocity, the mass loss during the SN explosion, and the current binary period. 
For HD 69686 with a long period (> 10 day s) and a very high ejection velocity (192 km s~^), 
the equation (7) in INelemans et al.l (119991 ) suggests a very large SN mass loss (> IOM0), 



which is large for a binary containing a late B-type star. Although the over-solar abundance 
of helium on HD 69686's surface could be a result of past mass transfer in a close binary. 



- 10 - 



it is not the only explanation because helium peculiar (both rich a nd poor) stars are quite 
common among young late B (single) stars (iHuang fc Gied l2006bl ). Also, as described in 
following paragraphs, we found evidence indicating that HD 69686 seems to belong to a 
co-moving cluster that naturally explains its high peculiar velocity. All of these arguments 
suggest that HD 69686 does not belong to a binary system, and we can estimate its age 
using single star evolutionary tracks. 

Muc h effort has been invested in searching for young massive stars forming in situ in 
the halo JZeliko fc Christodouloulll997l : llvnn et al.lbood iSooi : iMartinlbooi l2006h . Most of 
the blue halo objects examined in these studies are either ejected disk stars or old evolved 
stars (PHB or PAGE). Only a few stars are possible ca r ididat es of massive star formation 
in the halo. Ide la Fuente Marcos fc de la Fuente Marcos! (l2008l . and references therein) sum- 
marized the observational evidence of star formation far from the Galactic plane. One 
possible mechanism that can trigger star formation at such an unusual place is a col- 
lision bet ween high velocity c l ouds (HVC) and/or Galactic accretion fragments (or tidal 
streams). iDyson fc HartquistI (119831 ) proposed that early-type star formation can be trig- 
gered by "cloudlet-cloudlet collisions occurring every 10^ yr" in any particular high Galactic 
latitud e cloud, and the formation r ate could be 10^ early-type stars every 10* yr in the halo. 



Later, IChristodoulou et al.l (119971 ) suggested that the collisions between cloudlets within 
HVCs may occur less frequently than previously thought, and they derive a much lower star 
formation rate from the observatio nal data. More recently, another mechanisra of sta r for- 
mation in the halo was proposed by Ide la Fuente Marcos fc de la Fuente Marcos! (12008! ). i.e., 
that the tidal force of a passing massive cluster (such as a globular cluster) can trigger star 
formation in high Galactic latitude clouds. Though the proposed mechanisms are different, 
they all suggest a common scenario that a group of stars should be formed together in situ in 
the halo. If this indeed is the case for HD 69686, we expect to see a group of stars that share 
the same space velocity with and have the same age of HD 69686. Finding such co-moving 
stars would provide strong evidence to support their formation in situ in the halo. 

Because HD 69686 has traveled a very long distance in 73 Myr, it is expected that the 
co-forming stars, if they exist, would spread into a much larger space around the current 
position of HD 69686 due to the velocity dispersion in the original cluster. A dispersion of 
10 km s~^ in original space velocity would cause the co-forming group spread into a volume 
of 750^ pc^, i.e., spread all over the sky as seen from the Sun. To make the search viable, we 
limit our search range to the vicinity of the trajectory of HD 69686. Thus, our search will 
detect only a very small fraction of any co-moving stars whose space velocities are similar to 
that of HD 69686. 



To carry out this kind of search, we first need the projected trajectory of HD 69686 on 
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the sky as a baseline. Any co-moving objects are expected to be on or around this basehne 
and to have similar proper motions and radial velocities to the values of the nearest point 
on the trajectory. We calculated the projected positions of HD 69686's trajectory on the 
sky every 0.05 Myr between 0.5 Myr in the past and 0.5 Myr in the future. The results are 
given in Table 5 a nd plot t ed in Figure 5. We used the Catalogue of Stars with High-Proper 
Motions (V2)0 by llvanovl (120081 ) to search for similar high proper motion stars. The search 



region (about 2 degrees around the trajectory) is indicated as a shadowed area in Figure 5. 
We found a total of 26 stars (including HD 69686) from roughly 2000 high proper motion 
stars (> 40 mas yr^^) in the search region. Their proper motions (in both RA and Dec 
directions) have differences of 15 mas yr~^ or less from the reference points on the trajectory 
of HD 69686. All 25 stars other than HD 69686 are plotted in Figure 5 as open or filled 
circles with their proper motions indicated by line segments. 

Note that the proper motions illustrated in Figure 5 are not aligned with the projected 
trajectory (HD 69686 seems to move away from its own trajectory). The trajectory of 
HD 69686 in space (Fig. 4) is a curve in the Galactic potential field that differs from circular 
motion. In order to find the co-moving stars that move on a similar trajectory, we need to 
project the whole trajectory on the current sky. This means that both the past and future 
parts of the trajectory of HD 69686 have to be projected on the present sky. However, 
HD 69686's proper motion is related to projected positions on the sky at different moments, 
i.e., the past-time positions on the past-time sky and the future positions on the future sky. 
This is different from their position in the present sky (the solid line in Fig. 5) because 
of the Sun's orbital motion in the Galaxy. As the Sun (and Earth) is moving at ~220 
km s~^ around the Galactic center, the projected trajectory of HD 69686 (the solid line) 
will move downward (and slightly to the left) in Figure 5 (in the reverse direction of the 
solar circular motion, which is approximately parallel to Galactic latitude). Thus, although 
HD 69686 seems to move away from its projected trajectory on the sky, it actually stays on 
the trajectory because the projected trajectory is moving downward at the same time. 

Having similar proper motions does not necessarily imply that these stars belong to a 
co-moving group. Some foreground and background high proper motion stars can be mixed 
in. We can further sift these foreground/background stars out using the color-magnitude 
diagram (CMD) of the expected co-moving group. All members in the HD 69686 co-moving 
group are expected to have the same age (~73 Myr) and a distance around 380 pc (because 
we limit our search region in the vicinity of the trajectory whose distance is given in the last 
column of Table 5). A 70 Myr isochrone (for solar metallicity Z = 0.02 and log Age = 7.84) 



http://cdsarc.u-strasbg.fr/viz-bin/Cat7I/306A 
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extracted from evolutionary models by iLejeune fc Schaererl (120011 ) is plotted in Figure 6 
assuming a distance of 380 pc and Ay = 0. The JHK,^ magnitudes of all 26 stars were 
obtained from the 2MASS catal ogue (ISkrutskie et al.l |2006| ) and were transformed into the 
system of ISessell fc BrettI JlQSsh using formula^ (J - K)bb = [{J - ^)2M + 0.018]/0.983 
and i^BB = iKs)2m — 0.001 (J — K)bb + 0.039. From Figure 6, we can easily tell those 
foreground/background stars which fall far above or below the isochrone on the CMD. This 
step helped us to discard 13 out of 25 stars (the 13 deleted objects are plotted as open 
circles in Fig. 5 and Fig. 6). The remaining 12 stars plus HD 69686 are the final candidates 
of the high velocity co-moving group (plotted as filled circles in Fig. 5 and Fig. 6). The 
key parameters of the final candidates are listed in Table 6. However, we cannot be sure 
that all these 12 stars co-move with HD 69686 because we do not have the radial velocity 
data for them yet. New spectroscopic observations of these stars are planned in the near 
future. If some of the 12 stars are confirmed to co-move with HD 69686 by the follow up 
observation in the future, then we can rule out a runaway origin for HD 69686 because the 
mechanisms of generating runaways cannot account for the ejection of a group of stars in 
a common direction. Instead, the existence of a co-moving group would support the idea 
that the massive star formed in situ in the halo (most likely in a HVC or Galactic merger 
fragment) with other stars 73 Myr ago. 

The trajectory of HD 69686 covers a large range in all dimensions of the Galactic 
potential field (4 ~ 12 kpc in R, —1.8 ~ +0.9 kpc in z, and almost 180° of Galactic azimuth). 
Its shape may depend on the sel ection of different Galactic models. We did calculate the 
trajectory using different models (jPaczyhskilll990l : lAllen fc Santillanlll99ll ) and found similar 
results with trivial differences that are too small to change our analysis results. All the 
models used in our calculations assume a similar circular velocity a t a distance to Galactic 
center of 8 kp c (~ 2 20 km s~^). We note that some observations (IMivamoto fc Zhul Il998 



Mendez et al. 200G; Uemura et al. 2000; Shattow fc Loeb 2009; Reid et al. 



200% suggest 



that the circular motion of the solar neighborhood is faster than 220 km s~^ by 30~40 km 
s~^. Although it will be very interesting to see how much the trajectory of HD 69686 would 
change in a Galactic model characterized by a larger rotation speed at the solar circle, it is 
beyond the scope of the present paper. 
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Table 1. Radial Velocity Measurements 





K(Nl)^ 


K(N2)b 


Lines 


(km s"-*-) 


(km s-i) 


He I A4471 


149.0 


148.8 


Mg II A4481 


145.0 


148.2 


H I A4340 


149.2 


145.0 



^HJD(N1)=2,454,788.0079 
''HJD(N2)=2,454,789.0227 
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Table 2. Parameters of HD 69686 



Parameter 


Value 


/ (deg) 


214.5041 


b (deg) 


23.4482 


Vr (km s-i) 


148±6 


Vra (km s"^) 


-155±8 


Vr,ec (km s"^) 


13±1 


Vsini (km s^"*^) 


141±7 


Teff (K) 


14760±200 


logs' (dex) 


3.93±0.03 


log^poiar (dex) 


4.04±0.04 


Mass (Mq) 


4.42±0.15 


Age (Myr) 


73±10 


Distance (pc) 


380±20 


Ui^SK (km s-i) 


-191 


VisR'^ (km s~i) 


-37 


lyLSR" (km s-i) 


-57 



^The peculiar velocity of 
the Sun is Uq = 10.0, 
Vn = 5.25, and Wn = 7.17 
(iDehnen fc Binneylll998al ). 
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Table 3. Evolution Time of Rotating Stars'^ 



^/^crit 


Evolving Tinie*^ (Myr) 


0.1 


244 


0.3 


262 


0.5 


269 


0.7 


275 


0.8 


279 


0.9 


281 



'^This table is based on the data 
of t he rotating models o f 3 Mq 
from lEkstrom et al.l (120081 ). 

^The time that a rotating model 
takes to evolve from ZAMS to 

logS-polar = 4.04. 
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Table 4. Realistic Stellar Models at Different Inclination Angles'^ 



Inclination 


-^polar 


^polar 


log S-polar 


(degree) 


(K) 


(Re) 


(dex) 


90 


15260 


3.52 


3.99 


80 


15270 


3.52 


3.99 


70 


15310 


3.50 


4.00 


60 


15370 


3.48 


4.00 


50 


15440 


3.43 


4.01 


40 


15525 


3.32 


4.04 


30 


16260 


3.13 


4.09 


25 


16960 


2.95 


4.14 



^AU models have the same stellar mass 
(= 4.42M0) and same Vsini (= 141 
km s~^). The parameters of the models in 
this table are prepared in such way that, 
when we apply our H7 profile fit method to 
the synthesized H7 profiles of these mod- 
els, we obtain the same T^^r (= 14760K) 
and logfifmsr (= 3.93). 
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Table 5. Projected Trajectory of HD 69686 on the Sky 



Time 


RA 


Dec 


/ia cos S 




Vr 


Dist. 


(Myr) 


(J2000) 


(J2000) 


(mas/yr) 


(mas/yr) 


(km s~^) 


(pc) 


-0.50 


09 12 58.4 


-06 18 28 


-101.0 


26.7 


100 


389 


-0.45 


09 07 41.2 


-04 49 00 


-100.3 


25.6 


106 


386 


-0.40 


09 02 21.0 


-03 18 07 


-99.5 


24.3 


111 


384 


-0.35 


08 56 57.9 


-01 46 02 


-98.5 


22.7 


116 


381 


-0.30 


08 51 32.2 


-00 13 00 


-97.3 


21.0 


121 


380 


-0.25 


08 46 04.3 


+01 20 41 


-95.9 


19.1 


126 


378 


-0.20 


08 40 34.5 


+02 54 47 


-94.3 


17.0 


131 


378 


— U.io 


no Q c HQ O 
Uo OO 06. Z 


1 ri/l OQ ^7 
+U4 Zo 1 


— yz.o 


1/17 




Oi 1 


-0.10 


08 29 30.5 


+06 02 56 


-90.6 


12.3 


140 


378 


-0.05 


08 23 57.0 


+07 36 24 


-88.5 


9.8 


144 


378 


0.00 


08 18 23.0 


+09 09 05 


-86.2 


7.2 


148 


380 


0.05 


08 12 48.8 


+10 40 41 


-83.8 


4.6 


152 


381 


0.10 


08 07 14.9 


+ 12 10 59 


-81.3 


1.8 


155 


383 


0.15 


08 01 41.5 


+ 13 39 41 


-78.7 


-0.9 


158 


386 


0.20 


07 56 09.0 


+ 15 06 37 


-75.9 


-3.6 


161 


389 


0.25 


07 50 37.8 


+ 16 31 33 


-73.2 


-6.2 


164 


392 


0.30 


07 45 08.3 


+17 54 19 


-70.3 


-8.8 


166 


396 


0.35 


07 39 40.7 


+19 14 46 


-67.4 


-11.4 


169 


401 


0.40 


07 34 15.4 


+20 32 48 


-64.5 


-13.8 


171 


405 


0.45 


07 28 52.7 


+21 48 19 


-61.6 


-16.1 


172 


411 


0.50 


07 23 32.9 


+23 01 14 


-58.8 


-18.4 


174 


416 



Table 6. Candidates of the HD 69686 Co- Moving Group^ 



R A 






/\ / / r*r\G A 
LXjJjQ, LvUo U 








7 

o 








l^mas/ yr J 


i^mas/ yr J 


i^mas/yrj 


i^mas/yrj 


jTvei. 


t m Q rr 1 

l^magj 




08 18 22.98 


+09 09 05.1 


-86.2 


0.67 


7.2 


0.42 


1 


7.38 


7.51 


08 20 29.24 


+09 39 52.6 


-88.7 


5.9 


8.1 


5.6 


2 


13.60 


13.15 


08 19 28.74 


+08 26 23.8 


-79.4 




3.0 




3 


11.66 


11.37 


08 14 18.90 


+ 10 56 48.2 


-70.0 


36.0 


-6.6 


37.5 


4 


12.82 


12.41 


08 16 02.71 


+ 10 47 35.2 


-81.4 


28.9 


18.1 


30.5 


4 


12.68 


12.12 


08 05 53.09 


+ 13 58 28.0 


-93.2 




5.1 




3 


10.90 


10.83 


08 03 02.59 


+13 21 47.1 


-62.4 


5.5 


-13.2 


5.4 


2 


12.94 


12.58 


07 54 38.84 


+15 17 09.0 


-60.9 


2.1 


-14.9 


2.0 


2 


11.76 


11.44 


07 58 03.59 


+14 30 24.8 


-58.8 


5.2 


-8.9 


5.2 


4 


12.45 


11.78 


07 55 50.86 


+17 30 35.4 


-67.8 


5.8 


-19.7 


5.8 


4 


13.31 


12.47 


07 47 04.61 


+ 15 06 15.7 


-61.3 


6.0 


-11.2 


5.9 


4 


12.71 


11.91 


07 42 50.37 


+ 16 19 00.7 


-69.1 


6.1 


3.8 


6.0 


2 


13.01 


12.20 


07 58 14.04 


+ 17 18 01.8 


-61.5 


29.5 


-14.8 


31.0 


4 


11.78 


11.71 



''The first row is for HD 69686. 



^The Catalogue of Stars with High Proper Motions (llvanovi 120081 ) doe s not provide the e rrors 
in the proper motion s. T his column indicates sou rces for the error s : 1 - Ivan LeeuwenI (120071 ). 2 - 
Zacharias et~aD tooi ). 3 - iKislvuk et"aD J200oh . 4 - IPucourant et al.l J2006h . 
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Fig. 1. — The best fit solar abundance model (H:He = 0.9:1.0) of He I A4471 {solid line) 
compared to the observed spectrum {plus signs). The model parameters are: Teff = 16200 
K, logg = 3.92, V^sini = 141 km s'^. 




Fig. 2. — The best fit solar abundance model (H:He = 0.9:1.0) of the H7 profile {solid line) 
compared to the observed spectrum {plus signs). The model parameters are: T^s = 14500 
K, logg = 3.91, and Vsini — 141 km s~^. Note: the He I A4387 profile is not included in 
the synthesized spectrum. 
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Fig. 3. — The theoretical unreddened SED based on a Kurucz model ( ICastelli fc Kurucz 
20031 ) with Teff = 14760, logg = 3.93, and an angular diameter of 0.085±0.004 mas (limb 
darkened). The optical/IR (transformed from the B, V, Ic and JHKs magnitudes) and 
UV flux measurements are plotted as plus signs. 
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Fig. 4. — The trajec tory of HP 69686 (solid l ine) during the past 73 Myr based on the 
Galactic model 2 by iDehnen fc Binneyl (Il998bl ) (Kirc = 217 km s~^ at d = 8 kpc). The 
trajectories with slightly different initial velocities (by 1 a off in Vra) are plotted as dotted 
lines. The orbit of the Sun is plotted as a dashed line in the upper panel. The Galactic 
center is marked by the big plus sign. Diamonds indicate time intervals of 10 Myr. 
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Fig. 5. — The sky map of high proper motion stars around HD 69686. HD 69686 is plotted 
as a diamond. The candidate members of the HD 69686 co-moving group are plotted as 
filled circles. The foreground/background stars with similar proper motions are plotted as 
open circles. All oth er high prope r motion stars from the Catalogue of Stars with High- 
Proper Motions (V2) (Ilvanovll2008l ) are plotted as small dots. The solid line is the projected 
trajectory of HD 69686 on the sky. The plus signs mark the positions at 0.05 Myr intervals. 
The proper motion lines associated with individual objects indicate the size of their motions 
on the sky over a time period of 0.04 Myr. 
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Fig. 6. — The color-magnitude diagram of the HD 69686 co-moving group. HD 69686 
is plotted as a diamond. The candidate members of the HD 69686 co-moving group are 
plotted as filled circles. The foreground/background stars with similar proper motions are 
plotted as open circles. The solid l i ne is an isochrone of log Age = 7.84 and metallicity 
Z = 0.020 from iLejeune fc Schaererl (l200ll ) for a distance of 380 pc. The dotted lines are 
the same isochrone at distances 380±60 pc. 



